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EXISTENCE AND STABILITY OF STEADY-STATE SOLUTIONS 
TO THE QUASI-GEOSTROPHIC EQUATIONS IN 


MIMI DAI 


Abstract. We consider the stationary Quasi-Geostrophic equation in the 
whole space driven by a force /. Under certain smallness assumptions 
of /, we establish the existence of solutions with finite norm. This solution 
is unique among all solutions with finite energy. The unique solution © is 
also shown to be stable in the sense: any solution of the evolutionary Quasi- 
Geostrophic equation driven by / and starting with finite energy, will return 
to ©. 


1. Introduction 

In this paper we consider the existence and stability of stationary solutions to 
the two dimensional active scalar equation 


( 1 . 1 ) 


9t + u-V9 + kK°‘9 = /, 
u = R^9, 


in X (0,T) with T > 0, where a e [1,2), k > 0, A = V—A is the Zygmund 
operator, and 

R^e = K-^{-d29,di9). 

The scalar fnnction 9 represents the potential temperature and the vector function 
u represents the fluid velocity. The force / is independent of time. If 9t = 0, ini) 
reduces to a stationary equation: 


( 1 . 2 ) 


U-VQ + kA“0 = /, 

u = R-^e. 


Equation dm with a = 1 is a fundamental model of the surface quasi-geostrophic 
equation (SQG) [71[5S]. It describes the evolution of the surface temperature field 
in a rapidly rotating and stably stratified fluid with potential velocity. As pointed 
out in [7], this equation attracts interest of scientists and mathematicians due to 
two major reasons: it is a fundamental model for the actual geophysical flows with 
applications in atmosphere and oceanography study; from the mathematical point 
of view, the behavior of strongly nonlinear solutions to (ED with /t = 0, / = 0 
in 2D and the behavior of potentially singular solutions to the Euler’s equation in 
3D are strikingly analogous which has been justified both analytically and numeri¬ 
cally. For literature the readers are refereed to and the 

references therein. 

When a = 1, equation (ED is usually referred as critical SQG, although it is 
an open problem whether a dramatic change in the behavior of solutions occurs 
for the case of dissipation power less than I. The global regularity problem of the 
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critical SQG equation has been very challenging due to the balance of the nonlinear 
term and the dissipative term in (HU. This problem for the unforced critical SQG 
is resolved now by Kieslev, Nazarov and Volberg HI], Caffarelli and Vasseur [3], 
Kieslev and Nazarov [TQ] [20] and Constantin and Vicol [10] independently, using 
different sophisticated methods. 

The long time behaviors of the solutions to the critical SQG equation have been 
studied in [SlIHllgKIlKIlllMlEaET]. Estimates of the decay rates are obtained 
for mild solutions, regular solutions and weak solutions in Huns], HI], and [26] . 
respectively. In the forced case in |5] , with a special class of time independent force, 
the long time average behavior of viscosity solutions has been addressed and the 
absence of anomalous dissipation was obtained. In [18j . the authors proved that 
the linear instability implies the nonlinear instability in the energy norm in the 
periodic setting (on domain T^). Recently, Constantin, Tarfulea and Vicol [^ and 
Cheskidov and Dai |5] proved the existence of a compact global attractor in 
and L^(T^), respectively, with different assumptions on /. In particular, in [5], the 
assumption is that / is solely in for some p > 2. A crucial estimate is that 

any initial data Oq e L^(T^) yields weak (viscosity) solutions with bounded norm in 
L“, by using the De Giorgi iteration method. We point out that the boundedness 
in also plays an important rule in the stability analysis of the current paper. 

In this paper, we study the existence and uniqueness of solutions to (HU with 
finite energy {L^ norm) in the whole space and the nonlinear stability of such 
solution. On a domain where the Poincare’s type of inequality holds it is expected 
that the steady state has finite norm; and when the force is small, any evo¬ 
lutionary solution with finite norm converges to the unique steady state. But 
the Poincare’s type of inequality is not available in the whole space R^, where the 
situation is more delicate. Thus techniques other than the standard energy method 
are applied to study these problems. To establish the existence of solutions to HU 
with finite energy, inspired by the work of Bjorland and Schonbek [2] for Navier- 
Stokes equation, we make use of the fact: if $(t, x) is a fundamental solution for the 
operator dt + A°‘ = 0, then Jg is a fundamental solution for the operator 

A“ = 0. Therefore, the estimate of the steady state may be obtained through 
the fast decay estimates of solutions to the corresponding evolutionary equation. 
In the mean time, the Fourier splitting method (see HHl HI]) is at hand to estab¬ 
lish the decay estimates of the evolutionary solutions in the whole space R^. Once 
the finite energy estimate holds, the smallness assumption on the force yields the 
uniqueness of the steady state. To study the stability of the steady state, we com¬ 
bine the method of generalized energy estimates for low and high frequency parts 
of the solutions and the Fourier splitting method. 

In contrast to the NSE as studied in [2]) the SQG equation HU has a weaker 
dissipative term A“ with a e [1,2). This presents a serious obstacle to establish the 
stability of the steady state. However, thanks to the fact that a viscosity solution 
is bounded in L” (see |11|5]), this difficulty can be overcome. Thus, the stability is 
obtained in the sense: viscosity solutions of HU converge to the steady state 0. 

Before stating the main results, we hrst recall the definitions of weak solution 
and viscosity solution. 
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Definition 1.1. A weak solution to ()l.ll) is a function 6 e Cw([0, T]; 
satisfying that , for any cj) G x (0,T)), 

T T T T 

— r {9,4>t)dt— r {u94>) dt + K r {K'^9,K'^(j))dt = {9Q,4>{^^^))+ f {f,4>)dt. 

Jo Jo Jo Jo 

Definition 1.2. A weak solution 9(t) of (11.11) on [0, T] is called a viscosity solution 
if there exist sequences e„ —» 0 and 9n{t) satisfying 
f)G 

+ Un ■ V6(„ + KA°‘9n + CnAdn = f, 

R 9ji^ 

such that 9n9 vn Cw([0, T]; A^). 

Standard arguments imply that for any initial data 0o ^ there exists a viscosity 
solution 9{t) of (11.11) with a = 1 on [0,oo) with 9(0) = 9o (see [12], for example). 
The same result holds for the subcritical SQG equation with 1 < a < 2. It was 
proved in mis] that a viscosity solution to is bounded in A” provided f e 
for some p > ^. With this auxiliary estimate in A“, we are able to prove the 
nonlinear stability of the steady state 0. 

The main results for the critical SQG (a = 1) and the subcritical SQG with 
a G (1, 2) are stated separately as follows. 

Theorem 1.3. Suppose f e X = wi’'^ n A” satisfy the assumption: 

(A) f(f) = 0 for almost every |^| < pQ, for some constant po > 0. 

Then there exist a constant M > 0 and a constant C(po, n, M) so that if \\f\\x ^ 
C(po, K, M), system with a = 1 has a weak solution 0 G H'^, in the sense 

that for any </) e V, 

- (C/0, V</)) + K(Aiu,Ai&) = (/,</.). 

which satisfies 

(1.4) ||0||2<M, ||A^0||2 sg r-^II/IIx. 

Moreover, this solution is unique among all solutions satisfying (HD). 

Theorem 1.4. There exists a constant C(k) such that j/||/||jc ^ C(k), the steady 
state 0 obtained in Theorem \l.S\ is nonlinearly stable in the sense: let wq g A^ and 
9 be a viscosity solution of iny with 0=1 and initial data 9o = wq + Q which 
satisfies, for any T > 0, 

9 G A®(0, T; A^) n L^{0, T; H^) 

then 

lim \\9(t) — 0 II 2 = 0. 

>00 

Remark 1.5. The constant C(k) in Theorem 11.41 may be even smaller than the 
constant C(po, k, M) in Theorem 11.31 The smallness of ||/||a ^ CI(k) will imply 
the smallness of 0 in by Lemma 13.71 It is not known whether the smallness 
assumption on the force in Theorem 1 1.31 or Theorem [L4| can be removed or not. To 
apply the techniques presented in this paper, the smallness assumption is crucial. 
It is also worth to mention that after the work of [2] , the authors and their collab¬ 
orators proved the stability of the steady state of 3D NSE with rougher external 
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force in [1]. Namely, the assumption (A) was relaxed, but the smallness of / was 
still required. 

In the subcritical case a e (1,2), we establish that: 

Theorem 1.6. Let a G (1,2). Suppose f e X = n satisfying the 

assumption A. Then there exist a constant M > 0 and another constant C{po, k, M) 
so that */||/||x ^ C{po, K, M) the following hold: 

(i) System \1.2\) has a weak solution 0, in the sense that for any 4> B V, 

-{Ue,Vcf) + K{A^Q,A^f>) = {f,cf). 

And G e , with 

||0||2^M, ||At0||2=SK-i||/||x. 

Moreover, this solution is unigue among all solutions with finite norm in . 

(ii) The steady state 0 is nonlinearly stable in the sense: let wg G and 6 be a 
viscosity solution of the SQG equation G3P with initial data dg = Wq + G which 
satisfies, for any T > 0, 

Be L^{0,T;L^) n L'^{0,T] H'i) 

then 

lim \\9(t) — 0 II 2 = 0. 

»-00 

In these theorems and through the paper, we adopt the notations, || • ||p = || • ||lp, 
= Sr 2 fgdx and V = {(j) e |V • </> = 0}. 

We shall prove the results for the critical SQG (a = 1), that is, Theorem 1 1.31 and 
Theorem 11.41 For the subcritical SQG with a G (1,2), a similar analysis scheme 
with slight modification will yield the analogous result. Thus the proof of Theorem 
11.61 will be omitted. 

The paper is organized as follows: in section 2 we introduce the proof scheme and 
give some preliminary estimates; section 3 is devoted to the proof of the existence 
of steady state with finite energy, that is. Theorem II.31 and section 4 is devoted to 
the proof of the stability of the steady state, that is, Theorem ll.41 


2. Preliminaries 


2.1. Outline of the analysis scheme. The nontrivial part in the study of the 
stationary equation (HH) is to establish ||0||2 < M due to the fact that a Poincare 
type of inequality is not available in the whole space R^. Inspired by the work 
of [2] for the Navier-Stokes equation, we plan to achieve the goal in the following. 
Consider 


(2.5) 


— + U-V0 + KA°‘e = 0 

at 

0(0) = /, U = R^G. 


Formally, if 6 solves (12. 5L © = 9(t)dt solves 

G • V0 + kA“0 = / 
U = R-^G. 


( 2 . 6 ) 
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Since this equation is linear for a fixed U, solutions are unique and hence 0 = 0. 
Due to the integral Minkowski’s inequality, 

rCO rOO 

ii 0 ii 2 = ii wmhdt, 

Jo Jo 

one can see 0 e if || 0 (t )||2 ^ C'(l + t)“'*' with 7 > 1. Thus it is crucial to establish 
fast decay for 9 which will be addressed in Section [31 

2.2. Estimates on the operator A. We recall some standard estimates which 
will be used often through the paper. 

Lemma 2.1. Let ■^. If f satisfies Assumption (A), then 

||A‘'$(t)||p^e-"^o‘||A‘^/||p 

for any 1 ^ p ^ 00 and u ^ 0, where we adopt the convention hPv = v. 

For a proof of Lemma [331 see [3]. 

Lemma 2.2. |23j (Sobolev type inequality) Let 2 < p < co and u = 1 — There 
exists a constant C > 0 such that if v B S' and v is a function, then 

||z;||p ^ C||A"i;||2. 

Lemma 2.3. Let 1 < p < co. There exists a constant Cp depending only on p such 
that 

\\A''u{t)\\pikCp\\K'^e{t)\\p 

for all V ^ 0, t ^ 0 . 


3. Existence of Steady State with Finite Norm 

As outlined in Section|2l the fast decay of solution 9 to (12.51) as || 0 (t )||2 < C{1 + 
t)~^ with 7 > 1 plays an important rule. Due to a result in [11], the best decay 
one can expect for the critical case is that (since U does not decay, it is worse than 
the real SQG in decay rate) 

\\9{t)h^C{l + t)-^ 

through a standard Fourier splitting method (c.f. [551113]), where the initial data 
prevents a faster decay. To obtain a faster decay, we measure the difference 

/3 = 0- $ 

with $ = /, and /3 satisfies 

/3t + [/ • V/3 + kA“; 3 = -C/ . Vd> 

13 7) 

^ ^ /3(0) =0, U = E-^0. 

We expect that ||/ 3 (t )||2 decays as (1 + t)““ (see the proof for a = 1 in Section [3]). 
By assumption (A), the decay of $ is fast enough such that 9 will decay as fast as 

P- 

We first study the solutions for the two sequences of approximating equations: 

W ■ V0*+i + kA0*+i = / 

U' = R^0\ 


(3.8) 
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/3;+i + [/*. v/3*+i + kA/3*+i = -W ■ V$ 

(3.9) , ■ 

/3*+i(0) = 0, W = R^e\ 

For a fixed function W e we solve the two systems recursively to find ap¬ 
proximating solutions for (EH) and (Irfl) with a = 1. With uniform estimates on 
the approximating solutions, the convergence procedure will yield the estimate of 

||e||2 < M. 

3.1. Existence of Solutions to the Approximating Systems. 

Lemma 3.1. Let U^' e Lfi and f e X. There exists a unique weak solution 0®"''^ 
to \S. in the sense that for any </) G V, 

(3.10) (t/*0*+\ V<))) + k(A50*+i,A^) = 

Moreover, this solution satisfies 

(3.11) ||A50*+1||2^Ck-1||/||x. 

For the linear equation (13.81) , a standard Galerkin method gives the existence of 
solutions 0®+^ which satisfy 

k{ |A5 0®+i|2da;^ r /0®+ida; ^ C||A-5/||2||A50®+1||2 
Jk2 Jr2 

with an absolute constant C. 

Lemma 3.2. Let U'' £ with 

(3.12) \\Aiu%^CK-^\\f\\x 

and f satisfy assumption A. There exists a unique weak solution /3®+^ e L®(IR’'', L^)n 
L^(K.+ , H^) to \S. 9\) in the sense that for any (f) G C'^(R'''; V), 

-(/3®+\</)t) - (C/®/3®+\ V(^) + s:(A^/3®+i,A^) 

/3®+i(x,0) = 0, U^ = R^Q\ 

Moreover, it satisfies 

(3.13) sup ("||/3®+^(t)||^ + K r \\Ail3'^+'^{s)\\lds') ^CpQ^K~‘^\\f\\%. 

t>o V Jo / 

Proof: Applying Galerkin method yields a sequence of approximating smooth 
solutions which satisfy 

z at Jr 2 JJ 52 

= -(' (C/®-V$)/3)„+ida: 

Jr2 

<-||A^(C/®d>)||2 + ^||Ai/3)„+i||l. 

Ki Z 

If \Pl^^\^dx + K\ \Aipi;f^\^dx^-\\AHu^ml 

at Jr 2 Jr 2 k 


It follows 
(3.14) 
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Applying the Holder inequality and Lemma l2.11 we have 

||A5(C/*$)||1 ^ C\\K^Wm + C\\U^A^m 

^C||A5C/*||2||$||^ + C'||[/*||2||A^$||2 

^ C\\A^W\\l\\m, + C||A^[/*||^||A^$||^ 

Integrating over the time interval [0, t] for (I3.14p yields 

r \Pi^\t)\'^dx +K \ r \A^Pi^^\'^dx €CK-^\\f\\j^ \ e-^’^P°^ds 
Jr2 Jo Jr2 Jo 

< CPo^^-Wx- 

The inequality (13.131) follows by taking the limit to — > oo. 

□ 

3.2. Decay of /3L In this subsection we establish the decay for /3* by using the 
Fourier Splitting method. The rigorous proof can be carried through by working 
on a sequence of the Galerkin approximating solutions to (13.91) . To be abbreviate, 
we work on the regular solutions in a formal way. 

Lemma 3.3. Let be the solution of iS.tA) given by Lemma \3.‘A The following 
estimate holds 

(3.15) ll/3*+'(i)ll^<C'(po,«)||t/1l^||/||^(l + t)-' 
with the constant C{po,K,) depending on po and k. 

Proof: The first step is to establish 

(3.16) ^ C|^|||C/®||2 \\P^+%{s)ds + Po-iR-ill/b) . 

Indeed, taking Fourier transform of (13.91) yields 

/3®+i + • V^®+i) + ■ V$)) . 

Since /3®+^(0) = 0, we have 

/3®+i = - f • V/3®+i) + ■ V$)) {s)ds. 

Jo 

Due to the fact that V • ?7® = 0, it follows by the Young’s inequality and the 
Plancherel theorem 
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Then we claim that, for any I > 5, satisfies 

|[(i + tyii/3*+i||2] 

(3-17) ^c{i,po,^){i + ty-^\mi(^^\\p^+\s)hd^ 

+ C«-'ll/llA(ll<J’ll» + l|A^<i>0) (1 + t)', 


with the constant C{1, po, n) = C'k(^)®(1 + Pq for an absolute constant C. 

In order to prove it, multiplying (13.91) by and integrating over yields 


(3.18) -J/3^+^\\l + 4Aip 




A^C7*||2||$||^+ ||A5[/*||2||A2$||2 


< Cii ^WfWx (ll^ll® + l|A5d>||^ 


Denote S{t) as a sphere in with center at the origin and radius R{t) = -(l + t)“^. 
We decompose the viscous term in frequency space into two parts, corresponding to 
S{t) and its complementary S(tY, respectively. Applying the Plancherel theorem 
yields 



Combining ()3.18p and the last inequality it follows that 

^ |/3*+Tde + Cti-Wx (lei’ll™ + IIA^^IIl) . 

By (13.161) we have 


r i/3*+ipd^^c'iit/*i|2 

Js(t) 


ai+l 


s)\\2ds + PQ K II/II 2 


Sit) 






ai +1 


(S)||2ds+||/||2 R\ 


Thus, 

|||/3*+1||2 + nR\\/3^+m ^ CK\\W\\ji^ + p^^K-y \\P^+\s)hds + \\fh^ R^ 

Multiplying the last inequality by the time factor (1 + tY, it follows 

I [(1 + ty\\i 3 ^+Y\l] < cii,po, 8 :)iit/iii(i + ty-^ \w^^\s)\\2ds + 11/112) 

+ Cn-^WflYx (ll$l|» + I|A5$||2) (1 + ty 






STABILITY OF QUASI-GEOSTROPHIC EQUATION 


9 


with C{1, po, k) = C {1 + Pq^ for an absolute constant C, which concludes 

the argument of the claim. 

Combining (13.131) and (13.171) and using the fact ||/||2 < ||/||x yields 

^ c{i,po,>i){i + ty-nuy\i (tpo ^N-^ii/iii +11/112)' 

+ (||$||^ + ||A^$||2) (1 + ty 

^c{i,Po,^^)\mui + \\f\yx) wfVxii + ty-^ 

+ CK-Wx{\mi + \\A^m) {i + ty. 

Integrating in time and applying Lemma 12.11 it follows that 

{i + ty\\f5^+^mi 

^ C{1, Po, k)\\U^\\1 [1 + WfWx) WfWx f (1 + s)* ^ds 

Jo 

+ Ck ^II/II^J ^||<I’(s)||^ + ||A2(i)(s)|| 2 ^ (1 + s)*ds 

^ c{i, PO, n)\\u% (1 + ml) ml [( 1 + ty-^ - 1 ] 

+ cn-ymi (ii/ii» + w^hwi) £ + syds. 

One can verify that 

r e-2«Po*(l + s)'ds~ (l + t)'e-2'"^“* 

Jo 

which is bounded for all t > 0 since I, po, k > 0. It infers 

^ll/llx (ll/liro + 11^^ /II4) J e + ^)‘ds ^ C(l, Po, K)llflll. 

Thus, we obtain the preliminary decay estimate 

( 3 . 19 ) \\m\mi ^ c{i,po, mull (i + m\i) uwii^+mi 

Combining (13.171) and (I3.19L we repeat the procedure as above and obtain that 

\\m\t)\\mc{i,po,K)\\um2\\f\\i{i+t)-i 

□ 

3 . 3 . Proof of 0*+^ = m^dt solving (|3.8p . As pointed out in Section[21 for¬ 
mally 0®+^ = m^dt solves the approximating system (13.81) with 0®+^ = /3*+^-|-<I>. 
In this subsection, we prove this statement rigorously and show that 0®"''^ is uni¬ 
formly bounded in 

Lemma 3.4. Let he the solution of iS. .9|) given by Theorem 1 ,9. ill Then it 

satisfies e L^. As a eonsequence, is 

finite a.e. in K^. 
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Proof: For each i we define the sequence as 



We have 0^1^^ g since G for almost every t by Lemma [3.21 Applying 

Minkowski’s inequality, Lemma I^TT] and (13.151) yields 

r" ~ 

I|0n+'||2< 

JO 

nn nn 

^ \\i3^+\t)\\2dt+ wmhdt 

JO Jo 

^C{po,n)\\W\\l\\ffx+\\f\\x. 

Similarly, the Minkowski’s inequality implies 

rn + l 

Jn 

It is known that ||0®+^(t)||2(it is finite by Lemmas 12.II and 13.31 Thus, 

J 'Tl + l 

0 as n —> 00 . 

n 

We infer that the sequence {0^^^} is Cauchy in and hence it converges to a 
limit 9'^^^{t)dt e L^. 

□ 

Lemma 3.5. We have 0*+^ = that is, {t)dt solves \3.8\] with 

[/* satisfying \3.12\) and f satisfying the assumption A. 

Proof: Multiplying (13.9p by </> G V, integrating over space and time, we infer 
that from = 0*+^ — $ 

I + {AH^+\Aicl>)dt. 

Noticing that 0*+^(O) = /3*+^(0) — ^*(0) = /, exchanging the order of the integration 
yields 

(3.20) (0*+'(n),<^) + {AHw&+^),Ah) = -«(A^0),+\ A^) + (/,0). 

Since 6(*+^(n) = /3*+^(n) + $(n), it follows that the first term on the left hand side 
tends to 0 as n ^ 00 by (13.151) . On the other hand, denote (9*+^ = 
and we have 

\{Ai[U\e],+^ - 0*+i)], A^(/.)| ^ c\m\4&+^ - (=)*+i||2||V<^||4 

^C||A^17*||2||0),+^-0*+^||2||V</.||4 

|(Ay0),+i - 0*+i), a4)| ^ C\\Q4^ - 0*+1||2||A((.||2. 

Lemma[3)4]implies that the sequence {0)i’*'^} converges strongly to 0®+^ = 0'‘^^{t)dt 

in . Thus, we conclude from the last two inequalities that 
|(C/®-V(0),+i-(9®+i),<)>)|->0, n^cx) 

|(A^(0);^i-0®+i),A^)| ^0, n^oo. 
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Therefore taking the limit n ^ oo in (I3.20|) yields that 0®+^ is a weak solution of 
(13.81) . The uniqueness in Lemma [3.11 implies 0®+^ = 0®+^, which completes the 
proof of the lemma. 

□ 

3.4. Finishing the proof of Theorem 11.31 In this subsection we show that the 
approximating solutions 0* to (13.81) converge to a limit function 0 which solves the 
stationary equation (O with certain smallness assumptions on /. Additionally, 
we establish a uniqueness result in the sense stated in Theorem 11.31 

Lemma 3.6. Let 0*+^ be the solution of i3.8\) given by Theorem \3.1\ with W 
and f satisfying the assumptions of the theorem. Let M > 0 be a constant such 
that ||t /*||2 < M. Then there exists a constant C{po, k, M) such that if \\f\\x < 
C{po, K, M) then ||0*^^||2 < M. 

Proof: By Lemma 15751 Lemma [2.II and (13.151) we infer that 


(3.21) 

||0*+^||2 < 

PCX) 

Jo 

PCX) 

PCX) 



||/3*+^(t)||2dt + 

mm 



Jo 

Jo 



c{po.K)\m\mx 

+ WfWx 


Take M satisfying the quadratic equation 

C{po,k)M‘^Z‘^ + Z = M 

with a positive root Z = ^^ Af). It follows from (13.211) 

that if \\f\\x < C{po, then ||0*'''^||2 < M. By the boundedness of Riesz 

transform on it follows that ||17®+^||2 < M. 

□ 

The following lemma establishes the estimates for higher order derivatives of the 
sequence 0®+^. 

Lemma 3.7. Let 0*+^ be the solution of i3.3\} given by Lemma [3A\ with IT^ and f 
satisfying the assumptions of the lemma. Assume additionally that 0* e iJs with 

(3.22) ||A0*||2<Ck-1||/||x, ||Ai0*||2 < Ck-iI/Ix. 

Then there exists a constant C{k, M) such that if \\f\\x ^ C{k, M) then 0*+^ B Hi 
satisfying 

||A0*+1||2 ^ Ck-'II/IIx, ||Ai0*+i||2 < CK-^\\f\\x. 

The constant C{k,M) depends on k and M. 

Proof: Formally, multiplying ()3.8I) by A0®+^ (to make it rigorous we can use 
the Galerkin approximating solutions as test function here) yields 

k || A 0 *+ 1||2 = -{W ■ V0*+\ A0*+i) + (/, A0*+^) 

||C/*||ao||A0*+i||2 + CK-^\\f\\l + J||A0*+i||2. 
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Followed from the classical Gagliardo-Nirenberg inequalities [T^ by complex inter¬ 
polation (see also [22]), and by the boundedness of Riesz transform we have 

(3.23) ||C/*|U<C||Aic/*|||||C/*||| 

^C||Ai0*|||||0*|||. 

Combining the last two inequalities, the assumption p.22[) and Lemma 13.61 gives 
that 

||A0*+i||2 ^ C«-'||Ai0^|||||0*|||||A0*+i||^ + Cn-^f\\l 
^ CK-tM3||/|||||A0*+i||2 + C^-^fWl 

We choose C{k,M) such that ||/||a: < C{k,M) and Ck~^M^\ f\\]^ ^ 1/2. Thus, 

||A 0 *+i||i^C'/^- 2 ||/|| 2 . 

Similarly, we multiply (13.81) formally by A^0®+^ and infer that 
k||A50*+1||2 = _([/» . V0 *+\ A20*+1) + (/, A20*+1) 

^ |(A5[/* • V0*+\ A50*+1)| + \{W ■ A3V0*+\ A50*+1)| 

+ C«-i||A ^/||2 +J||Ai 0 '+i ||2 

^ Ck-^WA^U^ ■ V0*+1||1 + ||17*||oo||A50*+1||2 

+ C«-i||AVll^ + fl|Ai 0 *+i||^. 

By Lemma 12.21 and the boundedness of Riesz Transform it follows 
||A5C/*. ve^+^Wl ^ C\\A^U^\\l\\VQ^+^\\l 
^C||A 0 »|| 2 ||Ai 0 ^+i||l. 

Combining (13.231) and the last two inequalities yields 

||Ai0*+i||i ^ CK-^\\Aiff, + CAc-2||A0*||i||Ai0*+i||i 

+ CA-i||Ai 0 *||||| 0 *|||||Ai 0 *+i||l 

< CK-^Aif\\l + C («-"||/||i + «-tMi||/|||) ||Ai0*+i||i. 

Thus if we choose C{k, M) such that ||/||x ^ C{n, M) and 

c(«-4||/||i + «-tMi||/|||) ^ 1/2 

then 

\\AiQ^+^l^CK-Wx- 

□ 

We are now ready to prove the first main theorem. 

Proof of Theorem 11.31 Choose 0° G such that ||0°||2 ^ Af and ||A'''0°||2 ^ 
^~^II/I|a with 7 = 5,1,|. By the boundedness of Riesz transform on Sobolev 
spaces, we have 17° = i?-*-0° g 77°/^, ||C/°||2 ^ M and ||A'>'C/°||2 < «:~^||/||x with 
^ = 1 1 3 

Starting with 17°, we solve (13.81) recursively by Lemma [XT] and obtain a sequence 
{0*}”p with 17® = i?-*‘0®. It satishes that ||A'^ 0®||2 ^ k“^||/||x with 7 = 1, | 
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by Lemma [3.11 and Lemma [3.71 and ||0*||2 < M uniformly by Lemma [3.61 In the 
following we show that the sequence is Cauchy in and it converges 

to a limit function 0. 

Let = 0*+i — 0*. Thus solves the equation 

U^ . vr*+i + R-^Y^ ■ V0* + kAY^+^ = 0. 

Since Y^+^ e we can multiply the above equation by integrate by parts 

and we have, by noticing that V • C/* = V • R-^Y^ = 0 

K||A^y*+^||^ = -(V • {R-^Y^e^),Y^+^) 

^ c'K“^||A^(i?^r*0Oll2 + |l|A5y*+i||2. 

Thus by the boundedness of Riesz transform, Gagliardo-Nirenberg inequalities and 
Lemmas 13.61 and 13.71 we infer that 

(3.24) \\h^Y^+^\\l ^ C'K-2||A5(i?-Lr*0*)||2 

^ C'K“2||A5(i?-Ly*)0*||2 4 . C'K-2||i?-Lr*A5 0*||2 
^ C'k- 2||A5F*||2||0*|||^ + CK-2||y*||^||A^0*||^ 

^ C'K"2||A5F*||2||Ai0*|||||0*||| + CK“2||A5yi||2||A0*||2 

< + Cn-^f\\%) IIA^y^lli 

Applying the above estimate recursively we obtain 

||A^y*+i||2 ^ (ca-^||/|||m§ + CK-^f\\j,y WA^Ym 

^ {CK-f\\f\\lMi + C^-Wxy^-Wx- 

Therefore if we choose ||/||x small such that 

Cn-^\\f\\iMi+Cn-Wx<^ 

then T* ^ 0 in as i ^ 00 . It implies the sequence {0*} is Cauchy in and 
hence has a limit 0 e By Lemma [3761 we know ||0||2 < M. 

We briefly show that 0 is a solution to (II.2|) in the following. We take i ^ 00 in 
(j3.10|) . Since 0® —> 0 in iL 2 , it follows immediately that 

(A^0®+^,A^0) -> (A50,A^(/i). 

For the nonlinear term, we have 

{W-VQ^+\(l))-{U-WQ,(j)) 

= {{W - U) ■ V0®+\ (?!>) - ([/ • V(0®+^ - 0), (/)) 

= -((C/® - t/)0®+\ V(/)) + (C/(0®+i - 0), V(/)). 

By the boundedness of Riesz transform and Lemma 12.21 it follows 
|((t/®-C/)0®+\V</.)| ||C/®-C/||4||0*+'||4||V(^||2 

||A5(0®-0)||2||A50®+1||2||V</.||2. 

Since ||A50®||2 ^ Ck“^||/||x uniformly and ||V (/>||2 is bounded, we have 
((C/®-C/)0®+\V?i)0 asi^oo. 
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Similarly, we can show ([/(0*+^ — 0), Vcj)) 0 as * ^ oo. 

In the end we show that 0 is the unique solution of (jl.2|) among all the solutions 
satisfying 

(3.25) ||0||2<M, ||A50||2^Ck-1||/||x. 

Let 0 be another solution satisfying (13.251) . The difference T = 0 — 0 solves the 
equation 

(3.26) 17 • Vr + R-^Y ■ V0 + kAY = 0 

with U = R^Y. Since 0 and 0 are bounded in and Hi, we multiply (13.261) by 
Y to obtain that, by applying the procedure to derive (13.241) 

||A5y||2 ^ (cK-^WfWlAi^ + c^-^ffx) IIA^ril^. 

( 10 4 2 \ 

A CK~'^\\f\\x) < 1, which im¬ 
plies the solution is unique. It completes the proof of Theorem 11.31 


4. Stability of Solutions 


In this section we investigate the stability of the steady state solutions 0 obtained 
in Theorem 11.31 Specifically, we study the Quasi-Geostrophic equation (HID with 
the initial data Oq = 0 + wq, where wq is considered as a perturbation. It will 
be established that the solutions of (ED approach 0 as f ^ oo provided wq e L^. 
Namely, we shall prove Theorem [TH 

Let w = 9 — Q. By (11.11) and (11.21) w solves the equation 


(4.27) 


{wt + R^O ■ Vw + kAw + R-^w • V0 = 0, 
|ry(0) = Wq. 


A weak solution to (14.271) is a function w e C'w([0, T]; L^(IR^)) satisfying that , 
for any (jiG x (0,r)), 

1 1 

— {w,4't)dt + K {A^w,A^4’)dt— {R^9w,V(j))dt + {R^w ■ dt 

Jo Jo Jo Jo 

= (yjo,(f>(x,0)) + r {f,(j))dt. 

Jo 

We state the existence of solutions to (ED and (I4.27P as follows. 


Lemma 4.1. Let dg = 0 + wg with wg G L^(R^). Assume f e LP for some p > ^. 
There exists a viscosity solution 6 to im which satisfies, for any T > 0, 

9 G L®(0, T; ^ L‘^{Q, T;Hi). 

Moreover, 9 is bounded in ^“(O, T; L®(]R^)), that is, 

6(GL®(0,r;L“(R2)). 


Proof: Since 0 G L^(]R^) by Theorem 1 1.31 it follows e L^(]R^). The existence 
of viscosity solutions can be obtained by stand arguments (see [H]). The bounded¬ 
ness of 9 in L® on torus is obtained in [5l|4] (Lemma 2.3). An analogous proof 
will give the boundedness on (see [3]). 


□ 
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Lemma 4.2. Let 0 G There exists a constant C{k) such that if ||0||_f/i ^ 

C{k), the equation has a unique weak solution w satisfying the energy in¬ 

equality 


(4.28) 



< 




2 

2 - 


Proof: The existence can be obtained by the Galerkin approximating method, 
and the uniqueness follows from the linearity of the equation. To prove the energy 
inequality (14.281) , formally, we multiply (14.271) by w and use integration by parts to 
infer that 

TIT f \w\^dx-\-K{ \A^w\‘^dx = — f R^w-VQwdx. 

2 at Jr 2 Jr 2 Jr 2 

The right hand side can be estimated as 


R'^w ■ VQwdx 

Jr 2 


^ C'||i?^u;||4||ic||4||V0||2 


^ C||ia||2||V0||2 < q|A^u;||2||V0||2 


Thus if 0 G is small such that C'||V0||2 < n/ 2, we have 

{wl^dx + K^ ^ 0 

dt Jr2 Jr2 

which implies (14.281) by taking integration over the time interval (0, t) and taking 
the limit t —> oo. 

□ 


4.1. Proof of Theorem 11.41 We start with two generalized energy inequalities 
for the low and high frequency parts of w, respectively. 


Lemma 4.3. Let 4> = ,■)/)= 1 — , and function E{t) e C^([0,oo); L°^). 

Assume ||0||Lfi ^ C'(k). Then the solution w to Ii4-27\ l satisfies the following gen¬ 
eralized energy inequalities, 


(4 29 ) +2J {R -^0 

+ 2 r {R-^w ■ V0, ^ * ri;(s)) 

Js 


dr 

dr. 


(4.30) 


E{t)\\fjw{t)\\l !iE{s)\\'iljw{s)\\l - 2 k { E{t) |^|^u;(t) 

Js 

+ J E'{T)\\'iljw{r)\\2 dr2 ^ ^('^) 


dr 




R^0 ■ Vw, (1 - tP^)w 


R^w ■ V0, 
dr. 


dr 
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Proof: We only show the a priori estimates. Multiplying (I4.27|) by e 
(j) * w{s) and integrating over the whole space 

r * (j)* w(s) dx + k{ * cj) * w(s) dx+ 

Jk2 Jr2 

r R'‘~0 ■ * (/> * w(s) dx + f R'‘~w ■ ^ * w(s) dx = 0 

Jr2 Jr2 


The first two integrals can be written as 


1 d 

2 dt 


J 

Je 2 


^ — K(t — s)A 


^ * u;(s)p dx- ^ ^ wis)dt 

+K r * w{s)K 

Jr 2 ^ ^ 


* wdx 

* wdx 


l_d 

2 dt 


f 


^ — K{t — S 


* w{s) 


dx. 


Thus inequality (14.291) follows from the above equation by integrating over the time 
interval [s,f]. 

Taking Fourier Transform of (14.271) . multiplying it by ip'^wE(t), and integrating 
over yields 


1 d r 

2 dt Jr2 


E{t) \'il)w\^ “ F f d(, + KE{t) r 1^1 \iIjw\^ d^ 

^ Jr 2 Jr 2 

+E{t) I R^9 ■ VwiIj^w d^ + E{t) j R-^w ■ VQiIj^w d^ = 0 
Jr 2 Jr 2 


The inequality ()4.30p follows from the fact (^R^0 ■ VWjUi^ = 0 and integrating the 
above equation on time interval [s,t]. 

□ 

We use (j4.29|) and (j4.30|) to prove the low frequency and high frequency parts 
of w converge to 0 in respectively. The integrals on the right hand side of 
(I4.29l) - (|4.30l) will be estimated in the following. 

Recall that 0 is bounded in ^“((O, T) xR^), and || V0||2 (hence ||0||oo) is bounded. 
Thus w = 0 — Q is bounded in L®((0,T) x R^). Applying integration by parts, 
Holder’s inequality and Gagliardo-Nirenberg inequality, we have 


< 






{R^0 ■ Vrc, ^ * w'(s)) 

(j)* A“5 i^R^O ■ Vw) , w(s)^ 

(j)* (j) * A“2 (i?-*-0 • Vw) e“^”i'’’“®i^A2r(;(s) 

(j)* (j)* Ar^ {^R^0-Vw) A^w 
A'^ (^R^0') ||w||qo A^w + (j)*(j)*R^0 A^i 


A5, 






2 
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where we used the fact that, on the low frequency part, ^ ^ * R^d is bounded in 

L“°((0,T) X R^). Similarly, we have 


{R^w ■VQ,e * (j) * w{s)) 

{(j)* (j)* A“5cliv(i?-*-ii;0), 



A 2 div(i?"*-u>0) 

2 

A^w 


Ai{R^w)Q 

2 

A 

1 

+ 

2 


< 


K^{R^w) || 0 ||< 


^(||0||oo+||V0||2) 


A 2 


K^w 
2 
2 


i?-^wA^0 

+ ll^^wll 


A^ic 

A ^0 


A5^ 


w 


Therefore, it follows that 


limsup (j) * w{t) <(||0||oo + 
t^co 2 

+ limsup 
t^oo 


11V0112+C) r 

Js 


^-K.{t-s)A 


4> * w(s) 


A'^w 
2 

5 

2 


dr 


while it is known that limsupj_ 


e * w{s) 


= 0 since heat energy ap¬ 


proaches 0 as t —> 00 . By the energy inequality (I4.28p . the right hand side of the 
above inequality tends to 0 as s ^ 00 . Hence 


lim \\(j)w{t)\\l 

t^OO 


11 ^ ||2 

lim \\6 * I 

i —*•00 II II 2 


= 0. 


Next we estimate the contribution from the high frequency part. By the Holder’s 
inequality and Young’s inequality, we have 


R^wVe,ip^w^ = (a-i/ 2 V- 


< 


< 


A-1/2V • (i?-Lw0) 

2 


A^R-^ 


l|0IL 


=S(||0||«o+||V0||2) 


A^w 

A^w 


+ U? w 


A50 


A5, 


2 
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< 


< 


R^0 = (( 1 - 

(1 - |^|l/2 - 

2 

{l-ip^)J'iA^^^iR-^w)w) 

2 2 

+1 (1 - V’^)J■(Al/2(i?^0)^u)||^ I icr/^w)||^ 
+ (1 - 

2 


^||y^||oo 
=S(||V 0 b + C) 


A^w 

2 

+ 

2 

A^R^Q 

1^1 4 

4 

A^w 

+ c 

2 

A^w 


A 2 W 


where we used the fact that, the low frequency part (associated with the factor 
1 —'0^) of R-^9 is bounded in L”((0,T) x R^), since 9 is bounded in A”((0, T) x R^). 
Therefore, it follows from (I4.30|) that 


(4.31) 


'JW{t}\\l ||0ui(s)||2 


+ (||0||cx) + 


- 2 k r 

Js 

liveb + c) r 

J s 


m 

^ E{r) 

E{t) 


|C|2'0 w(t) 


dr + 


r* E' 

Js 'E 


E'(t) 9 

|| 0 u;(r )||2 dr 


A^w{t) 


dr. 


Let B{p) be a ball with center at origin and time dependent radius p(t) which will 
be determined later. Splitting the second and third terms in (j4.31() onto B{p) and 
B{py yields 


E{r) 

E{t) 


+ 


^ S 

Ejr) r 
Js E{t) J, 
r' E'jr) r 
Js E{t) J, 


|^|2 0 w ( t ) 


dr + 


f‘ E' 

Js 'e 


E'{t) 

E{t) 


= — 2k 

ri TA/ 


1^1 |?/;tD(r)| d^dr — 2K, 


ijw[r }\\2 dr 


f e{t) r 

Js E{t) J, 


B{pr 


|^||0u;(r)| d^dr 


< 


I 


B{p) 

* E'{t) — 2kE{t)p{t) 


|0'u;(r)| d^dT + 


f E'jr) f 
Js E{t) Jb(p)c 


|0ui(r)| d^dr 


E{t) 


r |0w(r)|^ + r \ |0'ui(r)|^ d^dr. 

Jb(p)<= Js 


Choose E{t) = (1 + 1)~^ and p{t) = 2 K{i+t) ^ certain large number 7 such that 

E'(t) — 2kE(t)p(t) = 0. 

Thus, it follows from (14.311) 




(4.32) 


(1 + sV 
(1 + t)^ 


w s 9 + 


I 


7(1 + r) 

(1 + t)' 


7-1 


f | 0 w(r)|" 

JB(p) 


+ 


+ IIV0II2 + C) A 2 u;(r) 


JB{p) 

dr. 


d^ dr 


2 
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Note '0(^) ^ for |^| ^ 1. So for large time 

r |'0 w(t)|^ dC < f |^|'^|w(t)|^ C(l + T)-'^||w(r)||2, 

JBip) Jb(p) 


>B(p) 

hence, for large s 

7(1 + r)^ 


f TTTT^ f |V'w(r)|^ d^dr ^ (7(1 + t) sup ||w(r)||2 

Js (1 + tp Jb(p) t^O 


>B(p) 

which tends to 0 as t —» 00 . Fix s first, then 


(1 + tp 


||i/''u;(s)|| 20 , as t —> 00 


since p is bounded and w 6 Therefore, (I4.32|) implies 

rCO 

limsup ||V''ui(t )||2 < (||0||oo + IV 0 II 2 + (7) A^i(;(r) 

t—>00 Js 

for arbitrary s. Taking s —>■ co yields 


dr 


lim WpwipW^ = 0. 

t—>co 


In the end, we conclude that 


lim ||w(t )||2 < lim ||((){(}(t)|L + lim ||i/;{(;(t)L = 0 

t—*<X) i—»-GO t—*^CC 


which completes the proof of Theorem 11.41 

Acknowledgment. The author would like to thank the anonymous referees for 
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